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The connectivity of the potential energy landscape in su- 
percooled atomic liquids is investigated through the calcu- 
lation of the instantaneous normal modes spectrum and a 
detailed analysis of the unstable directions in configuration 
space. We confirm the hypothesis that the mode-coupling 
critical temperature is the T at which the dynamics crosses 
over from free to activated exploration of configuration space. 
We also report the observed changes in the local connectiv- 
ity of configuration space sampled during aging, following a 
temperature jump from a liquid to a glassy state. 

PACS numbers: 64.70.Pf, 61.20.Ja, 61.20.Lc 



Understanding the microscopic mechanism for the in- 
credible slowing down of the dynamics in supercooled 
glass forming liquids is one of the hot topics in con- 
densed matter physics. In recent years, the combined 
effort of high level experimental techniques Q, compu- 
tational analysis [^|J|] and sophisticated theoretical ap- 
proaches has provided an enormous amount of novel 
information. In particular, it appears more and more 
clearly that in addition to the melting temperature and 
the calorimetric glass transition temperature, Tg, another 
temperature plays a relevant role. This temperature, lo- 
cated between 1.2 and 2 Tg depending on the fragility of 
the liquid, signals a definitive change in the microscopic 
processes leading to structural relaxation. Mode Cou- 
pling Theory (MCT) Q was first in identifying the role 
of this cross-over temperature, Tc. According to MCT, 
for T > Tc the molecular dynamics is controlled by the 
statistics of the orbits in phase space 0, while below Tc, 
the dynamics becomes controlled by phonon assisted pro- 
cesses Q . Studies based on disordered mean field p-spin 
models have also stressed the role of such a cross-over 
temperature [||. 

Computer simulation studies of realistic models of liq- 
uids have addressed the issue of the structure of con- 
figuration space in supercooled states Two differ- 
ent techniques have provided relevant information on 
phase space structure: the instantaneous normal mode 
approach (INM) ||], which focuses on the properties of 
the finite temperature Hessian, allowing the calculation 
of the curvature of the potential energy surface (PES) 
along 3N independent directions and (ii) the inherent 
structure (IS) approach which focuses on the lo- 



cal minima of the potential energy. Different from mean 
field models, computer simulation analysis provides a de- 
scription based on the system's potential energy, the free 
energy entering only via the equilibrium set of analyzed 
configurations. 

In principle, the INM approach should be well suited 
for detecting a change in the structure of the PES vis- 
ited above and below Tc. A plot of the T-dependence 
of the fraction of directions in configuration space with 
negative (unstable) curvature, /„, should reveal the pres- 
ence of Tc- Unfortunately, as it was soon found out [pr] , 
anharmonic effects play a non negligible role and several 
of the negative curvature directions are observed even in 
crystalline states, where diffusivity is negligible. A sim- 
ilar situation is seen in p-spin models, where the num- 
ber of negative eigenvalues of the Hessian is not zero at 
Tc Bembenek and Laird suggested inspecting 

the energy profile along the unstable directions to par- 
tition the negative eigenmodes in shoulder modes (sh) 
(i.e. anharmonic effects) and double-well (dw) modes 
(i.e. directions connecting different basins). While for a 
particular molecular system the fraction of double well 
modes fdw has been shown to go to zero close to Tc [0 , 
for atomic system, fdw is significantly different from zero 
even in crystalline and glassy states The existence 
of dw directions in thermodynamic conditions where the 
diffusivity is zero, i.e. in situations where the system 
is constrained in a well defined basin, strongly suggests 
that the two minima joined by the unstable double-well 
directions may lead to the same inherent structure IS. 
Such possibility has been demonstrated very clearly in 
Ref. [||. 

The aim of this Letter is to present a detailed eval- 
uation of: (i) the number of escape directions, Nescape, 
leading to a basin different from the starting one; (ii) the 
number of distinct basins, NdisUnct, which are connected 
on average to each configuration via a dw direction. This 
analysis, based on a computationally demanding proce- 
dure, shows that indeed the PES regions visited in ther- 
mal equilibrium above and below Tc are clearly different. 
We also study the evolution of the sampled PES as a 
function of time following a quench from above to below 
Tc- We show that, in analogy with the equilibrium case, 
two qualitatively different dynamical regimes exist dur- 
ing aging, related to different properties of the sampled 
PES. 

The system we study is composed of a binary (80:20) 
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mixture of TV = 1000 Lennard Jones atoms The 
dynamics of this system is weh described by MCT, with 
a critical temperature Tc equal to 0.435 JlTt - 

We begin by considering the equilibrium case. At each 
T, for each of the 50 analyzed configurations, we calculate 
the INM spectrum and — by rebuilding the potential en- 
ergy profile along straight paths following directions with 
negative curvature — we classify the unstable modes into 
dw and sh. It is important to notice that the classifi- 
cation is done by studying the shape of the PES along 
one eigenvector, i.e. beyond the point in configuration 
space where the eigenvector was calculated. In princi- 
ple, to identify a dw or a sh mode, one should follow a 
curvilinear path. Indeed the use of straight paths guar- 
antees only the identification of the dw modes whose one- 
dimensional saddle energy is close to the potential energy 
of the system. As discussed in more length in those 
dw modes are the ones relevant for describing motion in 
configuration space. 

Figure |] shows the T dependence of the unstable and 
dw modes for the studied system. In agreement with the 
analogous calculation of Ref. , even at the lowest tem- 
perature where equilibration is feasible within our com- 
puter facilities, fdw is significantly different from zero. 
To estimate roughly the number of double wells which 
do not contribute to diffusion, we have calculated the IS 
associated with the T=0.446 equilibrium configurations, 
and we have heated them back to various temperatures 
below Tc- The corresponding /„ and fdw are also shown 
in Fig.^ Both these quantities depend linearly on the 
temperature. If we extrapolate fdw for the equilibrium 
and for these non-diffusive cases, we find that the two 
curves cross close to the MCT critical temperature Tc- 
Thus Tc seems to be the temperature at which the num- 
ber of directions leading to different basins goes to zero, 
leaving only local activated processes as residual channels 
for structural relaxation |ll9|l . 

To estimate in a less ambiguous way the number of 
different basins which can be accessed from each con- 
figuration we apply the following procedure: (i) calcu- 
late the dw directions via INM calculation; (ii) follow 
each straight dw direction climbing over the potential 
energy barrier and down on the other side until a new 
minimum is found; (iii) perform a steepest descent path 
starting from the new minimum found along the dw di- 
rection, as first suggested by Gezelter et al. [|l^; (iv) save 
the resulting IS configuration; (v) repeat (ii-iv) for each 
dw mode. This procedure produces a list of ISs which 
can be reached from the initial configuration crossing a 
one-dimensional dw. We note that our procedure only 
guarantees that the two starting configurations for the 
quenches are on different sides of the double well, since 
— as discussed in — there is some arbitrariness in 
the location of the two minima. 

We calculate the relative distance da 
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between all IS pairs in the list to determine the num- 
ber of distinct basins, N distinct, connected to the starting 
configuration. Here a;^, y^, and are the coordinates 
of the m-th particle in the k-ih. IS of the list [ISk]- We 
also calculate the distance doi — where indicates the 
IS associated to the starting configuration — to enumer- 
ate the number of escape directions, Nescape, leading to 
a basin different from the original one (ISq). 

In Fig. 1^ (upper panel) we show a plot of the distribu- 
tions of the distances dgi between the starting IS and the 
minima identified with this procedure for different tem- 
peratures. At high temperatures, the distribution has a 
single peak centered approximately at d = 0.3. For lower 
temperatures, the peak moves to the left and decreases 
in height. At the same time, a second but distinct peak, 
centered around d w 10""*, appears. Since our sample 
is composed of 1000 particles, an average distance of the 
order of lO^'* means that, also in the case that only one 
particle has a different position between the two different 
configurations, this particle has moved less that 0.003 in- 
terparticle distance. Thus, we consider the IS i and j as 
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The upper panel of Fig.|| reports the T-dependence of 
the number of dw directions Ndw, the number of escape 
directions and the number of distinct basins which can 
be reached by following a dw direction. We find that at 
high temperature nearly all dw directions leads to a dif- 
ferent IS, thus fully contributing to the diffusion process. 
On lowering the temperature, a large fraction of the di- 
rections leads to the same minimum, and, close to 
almost all dw directions lead to intra-basin motion. A 
basin change becomes a rare event. Data in Fig.^ sup- 
port the view that Tc is associated to a change in the PES 
sampled by the system, and, consequently, in the type of 
dynamics that the system experiences. Above Tc the sys- 
tem can change basins in configuration space by moving 
freely along an accessible dw direction, while below Tc 
the thermally driven exploration of configuration space 
favors the exploration of the interior of the basins. Dif- 
fusion in configuration space, i.e. basin changes, requires 
an local activated process. 

We next turn to the study of out-of-equilibrium dy- 
namics. Configurations are equilibrated at high tem- 
perature (r = 5.0) and at time tw = are brought to 
low temperature Tf = 0.2 by instantaneously changing 
the control temperature of the thermostat. Within 100 
molecular dynamics steps, the kinetic energy of the sys- 
tem reaches the value corresponding to Tf. Configura- 
tions are recorded for different waiting times tw after the 
quench. We then repeat the same analysis done for the 
equilibrium configurations, as described before. 

Figure H (lower panel) shows the distribution of the dis- 
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tances doi between the minima ISq and I Si. For short 
the distribution has a single peak centered approximately 
at doi — 0.3, as found in the equilibrium configurations 
at high temperature. As the system ages, the peak shifts 
to smaller distances and a second distinct broad peak 
centered at approximately doi = 5 x 10^'* appears. The 
observed behavior is similar to the equilibrium one, on 
substituting T with i^. 

The lower panel of Fig. ^ shows Ndw, N escape and 
N distinct as functions of tu} ■ At small tw , following a dou- 
ble well direction always brings to a basin that is distinct 
from the original one, while at long times all dw direc- 
tions leads to intra-basin motion. Again, the behavior is 
very reminiscent of the equilibrium one, on substituting 
T with tw From these results, we suggest that two qual- 
itatively different dynamical processes control the aging 
processes, respectively for short and long t^^ For 
short waiting times the system is always close to saddle 
points that allow it to move from one basin to a distinct 
one. As the system ages, it is confined to wells that are 
deeper and deeper, and the number of basins that it can 
visit by following a dw direction decays to zero. Thus, 
for long the dynamics of aging proceeds through local 
activated processes, that allow the system to pass over 
potential energy barriers. As shown in Fig. |^, beyond 
tw = 10^ the number of distinct basins is almost zero, 
which, in analogy with the equilibrium results, supports 
the view that a cross-over from MCT-like dynamics to 
hopping dynamics may take place at a cross-over time 
during the aging process. 

In summary, this Letter confirms the hypothesis that 
Tc can be considered as the T at which dynamics crosses 
over from the free-exploration to the locally activated 
exploration of configuration space. In this respect, the 
analogy between the slowing down of dynamics of liquids 
and the slowing down of the dynamics in mean field Jo- 
spin model is strengthened. This Letter also shows that 
changes in dynamical processes associated with changes 
in the explored local connectivity of configuration space 
are observed even during the aging process, in the time 
window accessed by molecular dynamics experiments. 
MCT-based models of aging |^ could be able to de- 
scribe the early part of the aging process — i.e. the 
saddle-dominated dynamics — but may not be adequate 
for describing the locally activated dynamical region. 
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FIG. 1. T-dependence of /„ (squares) and fdw (circles). 
The filled symbols are calculated from equilibrium configura- 
tions, the empty symbols from out-of-equilibrium configura- 
tions below Tc (see text). 



FIG. 2. a) Distribution of the distances d = doi between 
the IS and the i-th IS from equilibrium configurations at 

different temperatures; b) same quantity for different t^^ val- 
ues during the aging process after a temperature jump from 
Ti = 5.0 to Tf = 0.2. 



FIG. 3. a) Number of double wells Ndw, of escape direc- 
tions Nescape and of distinct basins NdisUnct for equilibrium 
configurations as a function of T; b)Same quantities for dif- 
ferent tn, values after a temperature jump from Ti = 5.0 to 
Tf = 0.2. Time is measured in molecular dynamics steps. 
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